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ABSTRACT 


A  method  li  presented  for  predicting  the  orbit  of  a  passive  (non¬ 
transmitting)  satellite*  using  exclusively  Doppler  data  collected  during  one 
passage  of  the  aatellite  by  a  tracking  system  consisting  of  only  one  trans¬ 
mitter  and  one  receiver.  Equations  are  derived  which  relate  passive 
satellite  Doppler  tracking  data  to  range*  speed  and  time  parameters  at 
closest  approach.  These  determine  an  ellipsoid  whose  foci  are  the 
transmitter  and  receiver  stations.  The  satellite  lies  on  this  ellipsoid*  and 
ite  trajectory  is  tangent  to  the  surface.  The  satellite  position  and  velocity 
components  at  closest  approach  are  calculated  on  a  digital  computer  using 
information  contained  in  the  Doppler  data  together  with  the  geometry  of  the 
tracking  system.  Preliminary  orbits  determined  in  this  way  agree 
favorably  with  the  refined  orbits  of  BRL  and  Space  Track. 
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ORBIT  DETERMINATION  OF  A  NON -TRANSMIT TING  SATELLITE 
USING  DOPPLER  TRACKING  DATA 

1.  INTRODUCTION 

In  connection  with  the  application  q£  radio -Doppler  methods  to 
satellite  tracklngi  the  Ballistic  Research  Laboratories,  Aberdeen  Proving 
Ground,  requested  the  Space  Sciences  Laboratory  of  General  Electric  to 
Investigate  the  problem  of  orbit  determination  of  a  passive  (non -transmitting) 
earth  satellite  using  exclusively  Doppler  data  collected  during  one  pass  of 
the  vehicle. 

The  basic  problem  in  orbit  determination  is  to  determine  the 

components,  with  reepent  to  a  geocentric  inertial  coordinate  system,  of 

position  and  velocity  at  a  particular  time  for  the  orbit  being  observed, 

i.  e.,  determine  Jt(t  ),  y(t  ),  a(t  ),  x(t  ),  y{t  ),  a(t),  for  a  given  time 

o  o  o  o  o  o 

t^  .  For  an  earth  satellite,  this  is  equivalent  to  determining  the  orbital 

elements  a,  e,  i,  XI  ,  ,  M(t  ),  where 

o 

a  is  the  semi-major  axis, 

e  Is  the  eccentricity, 

i  is  the  inclination  of  the  orbital  plane, 

Di  is  the  right  ascension  of  the  ascending  node, 

00  is  the  argument  of  perigee, 

M(t^)  Is  the  mean  anomaly  at  t^  . 

The  equation  which  describes  the  Doppler  shift  of  radio  signals 
transmitted  from  a  station  on  the  earth,  reflected  from  the  passive 
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satellite,  and  received  at  another  station  on  the  earth  is* 

where  ^  and  jO  are  the  range  rates  of  the  satellite  relative  to 

the  transmitter  and  receiver  stations, 
respectively, 

"lip  and  are  the  transmitter  and  receiver  frequencies, 
respectively, 

and  c  is  the  speed  of  light. 


The  measured  Doppler  data  ("t  “  it  )  is  related  to  the  solution  x(t  ), 

T  K  Q 

^^*0^'  ^^*0^  of  the  orbit  determination  problem  by  the 
expressions  relating  the  time -dependent  range  rates  to  the  time -dependent 
solution  x(t),  y(t),  »(t),  x(t),  y(t),  2i(t)  of  the  equation  of  motion  of  the 
satellite,  A  "irMA/A  aC  ,  since  the  solution  of  the  equation  of  motion  is 


unique,  corresponding  to  a  particular  set  of  initial  conditions  x(t  ),  y(t  ) 

o  '  o 


z(to).  3l(t^). 

f^(t)=  (x(t)-] 

/(tg)*  The  desired  expression  for  ^  is 

«T.(t))  (x(t)  -X^(t))  +  (y(t)  -  Y^(t}  (y(t)  -  +  ( 

z(t)  -  Z^(t))  (z{t)  -  Z^(t)) 

(x(t)  -  Xp(t))  ^  +  1 

[y(t)-Yp(t))  ^  +  1 

^z(t)  -  Z^(t))  ^ 

where  X^,  Y^,  are  the  position  components  of  the  transmitter, 


Kp,  Y^,  are  the  velocity  components  of  the  transmitter. 

Similarly  for  • 

For  derivation  of  the  Doppler  equation  see  Appendix. 
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Having  obtained  Doppler  data  from  one  or  more  transmitter -receiver 
stations,  the  procedure  for  arriving  at  a  solution  to  the  orbit  determination 
problem  is  as  follows: 

1.  Select  preliminary  values  of  the  initial  conditions  y(*'Q)» 

z(t^),  establish  a  preliminary  orbit. 

2.  Solve  the  equation  of  motion  for  the  time -dependent  position 
and  velocity  components. 

3.  Compute  values  of  for  each  transmitter -receiver 

pair  at  times  corresponding  to  measured  Doppler  data. 

4.  Compare  computed  values  of  with  data. 

5.  Select  an  improved  set  of  initial  conditions  and  repeat 
the  process  until  the  computed  values  agree  with  the 
data. 

Items  2  to  5  constitute  the  determination  of  the  refined  or  final  orbit. 

The  Ballistic  Research  Laboratories  has  a  numerical  technique  for 
establishing  the  final  orbit  that  will  converge  with  data  from  a  single 
transmitter -receiver  system  provided  the  preliminary  values  of  the 
initial  conditions  are  sufficiently  accurate  [  1 

Doppler  tracking  of  active  (transmitting)  satellites  has  received 
considerable  attention  [2,  3,  4,  5,  6],  but  little  has  been  found  in  the 
literature  regarding  the  problem  of  establishing  preliminary  orbits. 

For  these  reasons,  emphasis  at  the  Space  Sciences  Laboratory  has 
been  placed  on  obtaining  a  reliable  method  for  determining  a  preliminary 
orbit  of  a  passive  satellite,  using  only  Doppler  data  collected  from  a 
single  transmitter -receiver  system  during  one  passage  of  the  satellite. 

The  method  is  formulated  in  this  report,  and  applied  to  actual  Doppler 
data. 

Numbers  in  brackets  refer  to  bibliography  at  end  of  paper. 
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The  geometry  of  such  a  system  can  cause  an  ambiguity  in  the  position 
and  velocity  components  of  the  satellite  at  a  given  time.  For  a  specific 
altitude,  the  satellite  can  be  at  one  of  four  possible  positions.  These 
positions  consist  of  two  pairs  of  points,  each  pair  consisting  of  points 
symmetric  with  respect  to  the  midpoint  of  the  base  line  between  the 
transmitter  and  receiver.  For  each  such  pair,  the  velocity  components 
are  uniquely  determined. 

It  is  shown,  however,  that  the  transmitter  antenna  beam  can  be 
oriented  so  that  ambiguity  in  the  direction  of  motion  of  the  satellite 
(signs  of  velocity  components)  can  be  completely  eliminated.  Further¬ 
more,  severe  restrictions  can  be  imposed  on  the  position  coordinates  by 
the  speed  of  the  satellite  and  the  location  of  the  satellite  relative  to  the 
antenna  beams  at  the  given  time. 

Thus,  ambiguities  caused  by  the  basic  geometry  can  be  considerably, 
if  not  completely,  eliminated  by  judicious  orientation  of  the  antenna  beams 
and  by  using  all  possible  information  contained  in  the  data. 

The  theory  is  formulated  analytically,  and  computation  can  be 
performed  on  a  digital  computer.  This  suggests  the  possibility  that  the 
orbit  of  a  passive  satellite  can  be  determined  using  only  Doppler  data 
from  a  single  transmitter-receiver  system  while  the  satellite  is  still 
visible  above  the  horizon. 

The  application  of  this  theory  to  active  satellites  is  also  described. 
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2.  ESTABLISHING  A  PRELIMINARY  ORBIT  FOR  PASSIVE  SATELLITES 


Integration  of  the  Doppler  equation  for  a  single  transmitter-receiver 


pair, 

==  At  d j  where 

'V  r  >  "  T  1?  > 

T 

(2.  1) 

yields 

r* 

(2.  2) 

which  indicates  that,  at  a  given  time  t,  the  satellite  lies  on  an  ellipsoid 
of  semi-major  axis  a(t),  with  foci  at  T  and  K,  the  sites  of  the  trans¬ 
mitter  and  receiver.  Since  T  and  R  are  known  locations,  the  eccentricity 
of  the  ellipsoid  is  obtained  from  the  expression  A.  0  »  .  Thus  the 

ellipsoid  is  compl'«tely  determined  once  the  constant  of  Integration 
(  Uq  known.  This  is  determined  in  the  following  manner: 

The  Doppler  equation  (2.  1)  may  be  written  in  vector  form  as 


ft  s, 


(2.  3) 


In  a  geocentric  Inertial  coordinate  aystom^neglectlng  the  earth's  rotation* 
V  ,  where  v  is  the  satellite  velocity,  and  equation  (2.  3)  becomes 


V«  f 

fT  fR 


(2.  4) 


Define  the  vector 


(2.  5) 


Then 

|pft)|  =  p(i)  -  (2.6) 

For  a  low  earth  satellite,  the  error  incurred  by  this  assumption  will  be 
small. 
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and  the  Doppler  equation  (2.  1)  becomes 

p(i)  ^ 

'  2 

or 

Elf  - 

P  '  2 


(2.  7) 

(2.  8) 


Define  the  time  of  closest  approach  as  the  time  at  which  is 

minimum  (  =  0  )  and  let  in  (2.  2)  be  this  time.  Then  the  constant 

of  integration  in  (2.  2)  is  the  major  axis  of  the  smallest  ellipsoid  upon 
which  the  satellite  will  lie*  the  ellipsoid  at  closest  approach.  From 
(2.  2)  and  (2.  6), 

i  TPr'*’  ^  ^  * 


Using  zero  subscripts  to  denote,  values  at  the  time  of  closest  approach, 
it  follows  from  (2,  4),  (2,  5)  and  (2.  8)  that 

y,  1  f. 

P«»  Pij  >  (2.  10) 

Further,  the  vector  p(i)  is  perpendicular  to  the  tangent  plane  of  the 
ellipsoid  at  the  location  of  the  satellite.  Hence,  both  vectors  ^  and 
lie  in  the  plane  which  is  tangent  to  the  ellipsoid  at  closest  approach,  at 
the  location  of  the  satellite. 

Thus,  the  following  conclusions  have  been  established: 

At  the  time  of  closest  approach  t^  the  satellite  lies  on  the  surface  of 
the  ellipsoid  of  semi -major  axis  ,  having  foci  at  T  and  R,  the  trans¬ 
mitter  and  receiver;  further,  the  satellite  trajectory  is  tangent  to  this 
ellipsoid.  The  location  of  the  satellite  on  this  ellipsoid,  together  with 
its  associated  velocity  and  the  time  t^  ,  will  determine  a  preliminary 
orbit  of  the  satellite. 
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Before  proceeding  to  determine  and  the  satellite  position  and 
velocity  components  at  t^  »  an  interpretatibn  of  is  in  order.  This 

vector  can  be  expressed  in  the  form 

where  a«  b,  c  are  time -dependent  scalars,  for  the  vectors  tT 

constitute  a  basis  in  three-dimensional  space.  The  vector  be 

interpreted  as  a  fictitious  velocity,  for  if  we  hold  the  base  of  pC'^) fixed  at 
the  origin  of  a  three-dimensional  vector  apace,  the  locus  of  the  end  point 
of  pet)  will  be  a  fictitious  trajectory.  We  take  as  origin  of  this  vector 
space  the  base  of  ,  i.  e. ,  p  at  closest  approach.  TMs  point  will  be  in 
the  plane  of  ^  and  on  the  normal  to  the  actual  trajectory  at  the  point 
of  the  trajectory  where  ^  minimxim,  and  at  a  distance  along  this 

normal  equal  to  1/2  this  minimum  vedue  of  vector  space 

there  will  be  a  sphere  of  radius  p,  =  ^  ^  ^mi'n  ficti¬ 

tious  trajectory  and  the  true  trajectory  are  both  tangent  to  this  sphere  at 
the  point  of  closest  approach.  Elsewhere  both  trajectories  are  exterior  to 
the  sphere. 


actual  trajectory 


r 

fictitious  trajectory 


If  we  now  assume  that  p  is  constant  (  =  )  in  the  neighborhood  of 

m  can  write,  in  this  neig 

pci)  s  +  i?o  (t-te)  ^2.  11)  * 


closest  approach  t^  ,  we  can  write,  in  this  neighborhood; 


Then,  using  (2.  10) 


and 


-  ’i 


(2.  12) 


(2.13) 
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Determination  of  Closest-Approach  Parameters.  -  If  the  measured  data 
consists  of  values  of  frequency  shift  at  specific  times,  equation  (2.  7) 
shows  that  the  parameters  t^,  p^,  can  be  computed  using  (2.  13).  If 

the  data  consists  of  average  values  of  frequency  shift  (  )  over  specific 

time  intervals  ,  the  mean  value  theorem  together  with  (2.  7)  and  (2.  13) 
gives 


A-f  /it  «  J 


(2.  14) 


and  the  desired  parameters  are  computed  from  this  expression. 

Great  care  must  be  exercised  in  fitting  the  right  side  of  (2.  13)  or 
(2.  14)  to  the  appropriate  data.  For  example,  (2.  13)  can  be  writte^n  as 


’  )/<t  rt  r» 


"i  “  i. 


(2.  15) 


rCT)  can  be  represented  by  the  series 


which  indicates  that  b/a  will  be  insignificant  if  the  data  used  to  deter"- 
mine  a  and  b  is  linear  inT.  This  can  lead  to  inaccuracies  in  the 
computation  for  b.  Similarly,  (2,  14)  can  be  written 

T,  -  4  ;  Tjj  s  t^-t^  ;  t,,  end  points  of  the  interval 


and  for  I  t-,  I,  /Tj  |  <  , 


*>  i,  Ccv  >• 
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Since  s  &  t  does  not  change  significantly,  it  is  apparent  that  if  the 

data  is  linear  in  the  discrete  variable  ,  then  again  b/a  will 

be  very  small. 

Since  the  use  of  (2.  13)  as  a  representation  of  the  Doppler  curve,  or 
S-curve  as  it  is  commonly  called,  is  restricted  to  the  neighborhood  of 
closest  approach,  this  equation  cannot  be  fitted  directly  to  the  measured 
data  if  noise  obscures  non-linearities  in  the  neighborhood  of  t^.  Further, 
formula  (2.  13)  is  symmetric  about  t^,  whereas  the  actual  S-curve  is 
asymmetric  for  large  values  of  t.  We  therefore  require  an  analytic 
representation  of  the  measured  data  which  is  non-linear  and  almost 
symmetric  for  small  T  =  t-t^,  to  which  the  symmetric  formula  (2.  13) 
can  be  equated  in  the  neighborhood  of  t^.  The  simplest  desired 
representation  is  a  polynomial  fit  to  all  the  measured  data. 

If  is  known,  standard  computational  procedures  can  be  applied  to 
fit  (2.  13)  or  (2.  15)  to  data.  Write  (2>  15)  as 


f.  =  , 


(2.  17) 


a  linear  expression  in  a  and  b,  and  the  usual  least-squares  formulas 
apply. 

Two  different  methods  have  been  used  successfully  to  determine  t 

o 

prior  to  fitting  (2.  15)  or  (2.  17)  to  the  analytical  representation  of  the 
data.  One  method,  easily  programmed  for  a  digital  computer,  is  to  fit 
the  cubic  in  t, 

n  s.  CL  i-  bt  f  j 

(2.  18) 

to  the  measured  data  in  the  neighborhood  of  t^,  with  the  time  origin 

for  t  fairly  close  to  t^  so  that  round-off  errors  will  not  be  significant. 

Once  the  coefficients  in  (2.  18)  are  determined,  t  is  computed  as  the 

o 

single  root  of  the  equation  p  =  o  that  lies  in  this  neighborhood. 
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An  alternate  procedure,'  which  has  been  very  successful,  but  more 

difficult  to  program  for  a  digital  computer,  is  based  on  formula  (2.  17), 

which  indicates  that  if  tV  p*  is  plotted  against  T*,  the  data  will  lie  on  a 

straight  line  with  slope  b  and  intercept  a.  The  parameter  t  must  be 

o 

known  in  order  to  measure  T.  of  course.  If  t  is  not  known,  denote 

o 

,  .i 

£  ^  %  '"fj,  >  t*  the  assumed  time  of  closest  approach,  and  let  be  the 
time  measured  from  ,  i.e.,  Tf  ® 4  ‘  Then  T  ■  Tf f  and  (2.  17) 
becomes 

X 

i> 


f  +  a 
/ 


f, 

a  » 


(2.19) 


Thus,  if  a  ^  is  selected  as  the  time  origin,  the  measured  values 
will  plot  against  as  a  straight  line  upon  which  is  super¬ 
imposed  the  curve  af;  the  straight  line  will  have  slope  b  and  intercept 
a,  as  before.  Clearly  f/m'fso,  which  suggests  that  different  values 
of  \  be  selected  until  one  is  reached  which  results  in  being 

linear  in  tf*  *  This  is  not  easily  adapted  to  digital  solution,  however. 
Instead,  we  proceed  as  follows. 

Denote  «  y  T*' 7  (2*  19)  becomes 


■f  =  ^  ‘ (aee  figure  1). 


(2.  20) 


Then 


and 


SQn  il  a  sgn  Cy^)  . 

^1* 


(2.  21) 


(2.  22) 


Now,  if  is  chosen  to  the  right  of  the  correct  value  t  ,  '-“t  >0  . 
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0<rl 


Selecting  values  of  ■?>*■’>  “ 'f/c  will  be  positive  and,  from  (2.  22), 

<0  .  Ifi  on  the  other  hand,  t.  is  selected  to  the  left  of  the 

correct  value  t^,  .  Again  using  values  oi  >  O  ,  /c  will  be 

negative  and  >0  •  Also,  it  is  seen  from  (2.21)  that  for  'P  ^  O 

and  finite,  ^  .-o  if  and  only  if  £  =  o  ;  further,  approaches 

zero  monotonically  from  above  (or  below)  as  converges  to  t^  from 

the  left  (or  right).  Similar  results  hold  if  values  of  are  chosen  in 

both  cases. 

The  procedure  then  is:  (1)  Choose  a  value  of  tj  ;  (2)  Select  three 
values  of  t  to  the  right  of  and  compute  the  second  difference  of 

at  the  center  value  of  t,  using  measured  p  data.  If  this 

.  i 

second  difference  is  negative,  say,  this  indicates  that  the  selected 

is  to  the  right  of  t^,  the  correct  value;  (3)  Choose  a  second  ,  to  the 

left  of  the  first  one;  (4)  Using,  for  convenience,  the  same  three  values  of 

t  previously  used,  compute  the  second  difference  of  y  again  at  the  center 

value  (the  values  of  are  now  different).  If  this  second  difference  is 

again  negative,  select  a  third  t,'  farther  to  the  left.  If  the  second 

,  < 

difference  is  positive,  select  the  third  slightly  to  the  right  of  the 
second.  In  either  case,  repeat  the  procedure  until  the  t/  is  found  that 
makes  the  second  difference  vanish.  By  selecting  the  three  values  of  t, 
used  in  computing  the  second  differences,  sufficiently  far  from  the  chosen 

.  I 

tf  ,  the  singular  situation is  easily  avoided. 

This  procedure  suggests  that  t^  be  determined  as  just  discussed,  and 
then,  using  (2.  17),  a  and  b  be  computed  as  the  intercept  and  slope, 
respectively,  of  the  straight  line  obtained  from  the  measured  data. 

Error  analysis  shows  that  the  parameter  a  can  be  found  in  this  way,  but 
not  b;  for  to  insure  accuracy  in  p^  to  one  mile,  and  in  (pd  to  100  ft/sec, 
it  is  necessary  that  time  measurements  be  accurate  to  .  4  seconds  in 
computing  the  intercept  a,  and  to  .00025  seconds  in  computing  the 
slope  b.  The  latter  requirement  is  clearly  too  restrictive  when  the 
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measured  data  consists  of  values  of  p  assigned  to  the  midpoints  of  .  15 
second  intervals. 

Instead,  the  most  successful  results  have  been  obtained  by  the 
following  procedure: 

(1)  Determine  t  ,  either  by  the  f -function  second 

o 

differences  (2.  14)  -  (2.  22)  using  measured  data, 
or  by  fitting  (least  squares)  the  cubic  equation 
(2.  18)  to  the  measured  p  data  in  the  neighbor¬ 
hood  of  t  . 

o 

(2)  Using  this  t  ,  use  least-squares  to  fit  the  cubic 

to  all  the  measured  data. 

(3)  Determine  a,  b,  and  p^  by  fitting  (lease  squares) 
equation  (2.  17)  to  values  of  p  given  by  (2.  23)  over 
a  sufficiently  large  neighborhood  of  t^.  The  inter¬ 
val  %  t  3C^of  the  time  interval  between  t^  and 
the  side  beams  has  been  found  adequate. 

Determination  of  Position  and  Velocity  at  Closest  Approach.  -  Once 
is  determined,  the  ellipsoid  at  the  time  of  closest  approach  is  known, 
since  the  eccentricity  is  ,  where  TR  is  the  length  of  the  base 

I 

line  between  the  transmitter  and  receiver,  the  foci  of  the  ellipsoid.  It 
remains  now  to  determine  the  location  of  the  satellite  on  this  ellipsoid, 
as  well  as  its  velocity  components. 

We  know  from  the  previous  theory  that  the  satellite  trajectory  is  tan¬ 
gent  to  this  ellipsoid.  Further  information  can  be  obtained  from  the 
Doppler  data  and  also  from  the  geometry  of  the  transmitter  and  receiver 

A  truncated  series  (2.  16)  or  a  truncated  Fourier  series  might  also  fit 
the  data  well. 
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beams.  In  particular,  the  DOPLOC*  transmitter  station,  located  at  Fort 
Sill,  Oklahoma,  consists  of  a  50  KW  transmitter  and  three  high-gain 
antennas.  All  of  the  antennas  can  be  excited  separately  from  the  trans¬ 
mitter  but  not  simultaneously.  Each  antenna  is  capable  of  radiating  a 
thick  fan-shaped  beam  having  dimensions  of  approximately  8  x  76  degrees. 
The  beam  from  the  center  antenna  is  directed  vertically  with  the  76 
degree  dimension  in  the  east-west  direction.  The  north  and  south 
antennas  are  identical  to  the  center  antenna  except  that  they  are  tilted  to 
radiate  fan-shaped  beams  20  degrees  above  the  north  and  south  horizons 
respectively.  This  arrangement  of  beams  is  intended  to  detect  a  satellite 
as  it  rises  over  the  horizon,  track  through  the  beam  width,  switch  to  the 
center  antenna,  track  through  its  beam  and  then  switch  to  the  third 
antenna  to  obtain  three  segments  of  data  from  the  characteristic  S-curve. 

The  DOPLOC  system  has  two  receiving  stations,  one  at  the  White  Sands 
Missile  Range  and  the  other  at  Forrest  City,  Arkansas. 

At  each  receiving  station  three  high-gain  antennas,  which  are  identical 
to  the  transmitting  antennas  except  in  power -handling  capability,  are 
tilted  to  see  the  same  volumetric  space  that  is  illuminated  by  the  trans¬ 
mitting  antennas.  Thus  each  center  receiving  antenna  has  its  fan-shaped 
beam  in  the  east-west  vertical  plane  but  tilted  toward  F'ort  Sill.  In 
addition  to  being  tilted,  the  north  and  south  receiving  antennas  at  each 
receiving  station  are  rotated  so  that  each  station  sees  approximately 
half  of  the  volume  of  space  illuminated  by  the  north  or  south  transmitting 
antenna.  The  receiving  antennas  are  switched  synchronously  with  the 
transmitting  antennas. 

In  the  following  analysis,  only  one  transmitter -receiver  pair  will  be 
considered.  Without  loss  of  generality,  we  will  describe  the  technique 
for  the  transmitter  at  Fort  Sill,  Oklahoma,  and  the  receiver  at  Forrest 

*The  word  "DOPLOC”  is  a  combination  of  the  first  three  letters  of 
Doppler  and  lock  and  implies  a  phase -locked  Doppler  measuring  system. 
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City,  Arkansas.  These  stations  are  separated  by  a  base  line  of  435 
statute  miles,  and  of  course  are  the  foci  of  the  ellipsoid  of  closest 
approach  for  each  satellite  observed.  We  will  assume  that  this  ellipsoid 
(  Pq  )  has  already  been  determined. 

A  typical  Doppler  data  plot  for  this  transmitter -receiver  pair  is  shown 

in  figure  2.  Here  the  time  of  closest  approach  occurs  while  the  satellite 

/ 

is  in  the  center  beam,  which  means  that  the  satellite  lies  on  the  associ¬ 
ated  ellipsoid  of  closest  approach  in  the  strip  intersected  by  the  center 
beam  of  the  transmitter.  In  all  cases,  the  ellipsoid  of  closest  approach 
associated  with  a  particular  transmitter-receiver  pair  will  be  intersected 
by  the  three  transmitter  beams,  and  the  location  of  the  satellite,  at  time 
of  closest  approach,  relative  to  these  beams  gives  a  good  indication  of 
where  the  satellite  lies  on  the  ellipsoid.  Further,  the  time  of  appearance 
in  the  beams  indicates  the  general  direction  of  the  satellite.  For  example, 
the  satellite  of  figure  Z  is  clearly  going  south  as  it  passes  over  the 
stations. 

This  information  can  be  expressed  analytically,  in  a  form  easily 
programmed  for  a  digital  computer,  in  order  to  calculate  the  approxi¬ 
mate  position  and  velocity  components  of  the  satellite  at  the  time  of 
closest  approach. 

At  the  time  of  closest  approach,  a  horizontal  plane  at  the  altitude  h 
of  the  satellite  above  the  transmitter  will  intersect  the  transmitter  beams 
as  shown  in  figure  3.  If  we  assume  that  the  satellite  remains  in  this 
plane  during  its  flight  through  the  three  beams,  and  also  that  its  speed 
or  is  that  corresponding  to  a  circular  orbit  at  this  altitude  h,  then  the 
satellite  path  length 4S  between  the  center  and  north  beams  is  given  by 

AS  =  V  At  (2,  24) 

where  At  is  the  time  interval  corresponding  to  the  gap  in  the  data 
between  the  center  and  north  beams.  The  acute  angle  G  between  the 
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FIGURE  2- TYPICAL  DOPLOC  DATA  PLOT 
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FIG.3-  ORIENTATION  OF  ELUPSOID  AND  ANTENNA  BEAMS 


path  and  the  vertical  plane  of  the  transnnitter  center  beam  (in  this  case 
this  plane  is  in  the  east-west  direction)  is  obtained  from 


- 


0<  9^ 


s/n9  '  -  (2.25) 

where  juN  is  the  perpendicular  distance  between  these  beams  at  the 
given  altitude  h. 

It  will  be  convenient  to  use  a  topocentric  coordinate  system,  having 
origin  at  T,  x-axis  vertical,  z-axis  in  the  horizontal  plane  at  the 
elevation  of  the  transmitter  and  in  the  direction  of  the  ellipsoid  base  line 
TR  extended  westward;  y-axis  in  the  same  horizontal  plane  and  normal 
to  TR  (figure  3). 

Then,  corresponding  to  a  given  altitude  h  (or  x),  equations  (2.  24) 
and  (2.  25)  give 

0 , 


Sin 


fit 


(2.  26) 


where 


A 

r 


cot  28°  -  tan  4° 


-  radius  of  earth  at  transmitter 
k  =  Newton's  gravitational  constant 

M  =  mass  of  earth. 

Clearly  the  altitude  x  of  the  satellite  will  be  restricted  by  the  upper 
bound  of  sin  0  .  This  restriction  can  be  seen  intuitively:  the  higher  the 

satellite,  the  greater  must  be  its  speed  to  traverse  the  wider  distance 
between  beams  in  the  fixed  time  interval;  on  the  other  hand,  the  speed  in 
a  circular  orbit  reduces  with  increased  altitude.  There  will  thus  be  an 
upper  bound  on  permissible  altitudes. 

Now  denote  by^u  ,  0  i ytt  <  S  ,  the  angle  at  T  measured  clock¬ 

wise  in  a  horizontal  plane  from  the  base  line  TR  to  the  vertical  plane  of 
the  transmitter  center  beam;  and  by  o<  ,  —yx.  <  «  <r  TT-yu,  the  angle  at 
T  measured  counterclockwise  in  a  horizontal  plane  from  the  base  line 
TR  to  a  line  through  T  parallel  to  the  satellite  path  (figure  3). 
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In  the  horizontal  plane  at  some  altitude  x,  the  satellite  will  lie  on  the 
ellipse 

4  JL  =/ 

(2.27) 


where  a,  b  are  the  semi-major  axis  p^^  and  semi-minor  axis  respec¬ 
tively  of  the  ellipsoid  at  closest  approach.  Since  the  satellite  path  is 
tangent  to  this  ellipse  at  the  location  of  the  satellite,  this  point  is  to  be 
determined  from  the  expressions 

1  a 


i.’  7-  i««'« 
a* 

J  a'--*’)  _ 

k  f  tcin*(X 


(2.  28) 


(2.  29) 


The  velocity  components  of  the  satellite  are  to  be  determined  from  the 
expressions 


^  -0 

SIS  ^ 

2*  s  a'*  C6<  (X 


r+-  X 


(2.30) 


Clearly,  the  problem  is  not  yet  solved,  first  because  is  not  yet 
known,  and,  in  addition,  we  must  establish  a  criterion  for  determining 
the  signs  of  y,  z  +  ae,  y  and  z. 

From  figure  3  it  is  clear  that  for 


yu.<  0<S  y,  , 

(2.31) 

0(  is  given  by 

(2.32) 

For 

IL -yu, ^ 

(2.33) 
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(2.  34) 


Ot  is  given  by  ~  Tr~~  (O  iyH’)  . 

Due  to  the  manner  in  which  0  is  determined  (2.  25),  we  do  not  know 
which  expression  is  the  correct  one  for  Of  .  We  do  know,  however,  that 
all  possible  paths  are  uniquely  defined  either  by  or  by 
From  (2.  32)  and  (2.  34), 


if 

C  <  6 

then  -yu.  <  ^  0  and 

n-  ^ 

if 

yU<  9c  y 

y- 

then  iZ-2/a  and 

r  2  X 

if 

JL 

2 

then  i  E-^  and 

JLyc  1  . 

Further,  it  is  clear  from  figure  3  that 

if  ®<<  0  ,  the  satellite  lies  in  the  1st  quadrant  (Qj^)  or  3rd 

quadrant  (Q^)  of  the  ellipse  (2.  27); 

if  0<  <X<  ^  ,  the  satellite  lies  in  the  2nd  or  4th  quadrant  (Q2  or  Q4) 

if  ,  the  satellite  lies  in  the  Ist  or  3rd  quadrant  (Qj  or  Q3). 

The  relationships  that  exist  between  the  computed  values  of  ^  ^  , 

and  ofjj.  and  the  quadrants  in  which  the  satellite  might  lie  are  shown  in 
figure  4,  In  general,  a  given  value  of  &  determines  four  possible 
points  on  the  ellipse,  consisting  of  two- pairs  of  symmetric  pointSi  A 
unique  velocity  is  associated  with  each  pair.  If  all  four  points  lie  in  the 
1st  and  3rd  quadrants,  then  the  direction  of  the  satellite  is  clear,  i.  e. , 
y/z  >0  .  Further,  the  sign  of  y  is  determined  by  observing,  from  the 
data,  which  beam  the  satellite  crosses  first.**  Hence  the  sign  of  z  is 
known.  Similarly,  if  all  four  points  lie  in  the  2nd  and  4th  quadrants, 
y/z  <0  . 

If  the  center  beam  is  rotated  above  the  base  line,  and  if  yu  is  measured 
counterclockwise  from  the  base  line  to  the  beam,  the  formulas  will  be 
different  but  the  results  will  be  the  same. 

The  satellite  of  figure  2  is  travelling  in  a  southerly  direction,  hence 
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FIG.4-RELATI0NS  BETWEEN  SATELLITE  DIRECTION  AND  ITS  LOCATION 


Thus,  the  following  conclusions  have  been  established: 

For  a  given^  ,  and  altitude  x,  and  hence  0  (2.  26),  and  denoting 


B.  -  ^ 


^  +  iary*(By^\ 
La^ 


b-A 


tan(9yi)’ 


c  ■=  k 

cc 


r  «  4  ^ 

^  i'  irein*‘(&y^) 

L  ^  } 

tain(&y^  } 


U  0<  &  $ »  there  are  two  possible  points  in  1st  quadrant  and  two 
<C'  C.  C*)  possible  points  in  3rd  quadrant.  Since  V/z  >  0  for 
^  ^  all  these  points,  the  general  direction  of  mption  of  the 

satellite  is  known. 

Points  in  Ist  quad:  y  =  Bp  z  =  -ae  +  B2;  y  =  T  V  sin  (  Qyu)*  i  =  t  f  cos 

y  =  Cp  z  =  -ae  +  Cg)  y  =  ±  X''  sin  ( 6?^),  z  -  ±  x’’  cos  {Q+yi, 

Points  H  '^rd  quad:  y  =  -Bp  z  =  -ae  -B2:  y  =  T  x*"  sin  (  0yu. ),  z  ~  ±  V  cos  (  0yu. 

y  =  -Cp  z  =  -ae  -C2;  y  =  ^  X'‘8in(  ),  z  =  is  X*"  cos  (  Qy*- 

If  yt**  <  '  there  is  one  possible  point  in  each  of  the  four  quadrants. 

C  ^  direction  of  motion  of  the  satellite  is  not  known.  ** 

*The  sign  of  y  is  to  be  determined  frorri  the  data  for  each  satellite. 

**Although  the  sign  of  y  is  known  from  the  data,  the  sign  of  z  depends 
upon  whether  the  satellite  is  either  in  the  1st  or  3rd  quadrants,  or  in 
the  2nd  or  4th  quadrants.  This  is  not  known  in  this  case. 
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Point  in  1st  quad: 
Point  in  2nd  quad: 
Point  in  3rd  quad: 
Point  in  4th  quad: 


y  =  Cj^,  z  =  -ae  +  C2:  y  =  i  V  sin  ( 0  iyi),  z  =  d:  v'  cos  ( 0+^ ) 
y  =  z  =  -ae  +  B2;  y  =  ±  -v'  sin  (  z  =  ^  V  cos  {6pu.) 

y  =  -C^,  z  =  -ae  -  Cg,  y  =  ±  a'"  8in(  Q-yi),  z  =  ±  ir  cos  {By*.) 
y  =  -Bj,  z  =  -ae  -  B^:  y  =  *  Vsin(@^),  z  =  T  V  cos 


If 

^  /  2 
(B^>o,  C,>o) 


there  are  two  possible  points  in  the  second  quadrant  and 
two  possible  points  in  the  4th  quadrant.  Since  ^/z 
for  all  these  points,  the  general  direction  of  motion  of 
the  satellite  is  known. 


Points  in  2nd  quad:  y  =  Bp  z  =  -ae/+  B^:  y  =  *  V  sin  (0^ ),  z  =  =F  I'^cos  (  £?^) 

y  =  Cp  z  =  -ae  +  C2!  y  =  ±  V  sin  {Qpu-),  z  =  ±  V  cos  (  Oiyu.) 

Points  in  4th  quad:  y  =  -Bp  z  =  -ae  -  B2i  y  =  ^  "Lf  sin  ( 0^),  z=TVcos(0^) 

y  =  -Cp  z  =  -ae  -  C2J  y  =  dr  •V8in(0'^),  z  =  ±  Trco8(0^), 


Clearly,  the  ambiguity  of  direction  in  the  second  case  (  yj,  <  B  <  ^ yi  ) 
is  removed  if  yj,  =  ^  ,  i.  e. ,  if  the  transmitter  antenna  is  rotated  45°  from 

the  base  line. 

tt 

We  still  have  the  problem  of  selecting  the  correct  altitude  x,  and  the 
correct  point  (  y,  z  )  at  this  altitude,  out  of  four  possible  values  of  (  y,  z  ). 
We  have  seen  that  the  possible  valuer/  of  x  are  already  restricted  by  the 
upper  bound  of  sin  Q  (2.  26).  For  a  chosen  x,  the  four  possible  values  of 
(  y,  T,  )  just  defined  are  not  all  acceptable.  The  data  tells  us  where  the 
satellite  is,  at  t^,  relative  to  the  center  beam.  That  is,  the  satellite  is 
either  in  the  beam  or  on  one  side  of  it,  and  proceeding  toward  it  or  away 
from  it.  This  information,  expressed  analytically,  imposes  severe 
restrictions  on  the  possible  values  of  (  y,  z  )  that  are  to  be  accepted. 

We  first  construct  a  new  coordinate  system  x',  y\  z' ,  by  rotating 
the  y  and  z  axes  of  figure  3  clockwise  about  the  x-axis  through  the 
angle  M.  .  The  x'y'  plane  is  then  vertical  and  normal  to  the  center  beam 
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(figure  5).  Denote  by  tj^  and  the  times,  at  which  the  satellite  crosses 
the  edges  of  the  center  beam.  The  distaiice  travelled  by  the  satellite  along 
its  path  in  the  horizontal  plane  at  altitude  x  =  x'  between  the  time  of 
closest  approach  t^  and  the  time  it  crosses  the  vertical  plane  x'z'  Is 
given  by 


V 


- 1 


Then,  for  the  case  where  ^'>0  (satellite  travelling  north),  the  y'  co¬ 


ordinate  of  the  satellite  at  t^  is  given  by 


y;  =  V  s«  e(t.  -  ^ 

I  ^ 

and  for  ^  <0  (satellite  moving  south),  the  y  coordinate  at  t^  is  given 


Since  ~  ^  M  siy4 ,  the  appropriate  formula  can  be  used  as  a 

criterion  for  selection  of  permissible  points  (  y,  z  )  at  a  given  altitude,  x, 
from  the  four  points  previously  computed.  To  allow  for  possible  inaccuracy 
in  the  precise  location  of  the  satellite  at  times  t^,  t^i  and  t^ ,  we  can 
accept  values  of  the  y'  coordinate  that  lie  in  the  range  f  <  i  >0* 

To  Select  this  range,  define  the  angle  (f)  in  the  x'y'  plane  by  the  relation 


Then,  assigning  the  small  angle  A(p>0  ,  the  x,  y,  z  coordinates  of  the 

satellite  must  satisfy  the  inequality 


(2.  35) 


Application  to  Doppler  Field  Data.  -  Using  the  methods  just  described, 
preliminary  orbits  were  established  for  two  revolutions  of  the  Discoverer 
XI  satellite.  The  data  was  obtained  from  the  ARPA  -  BRL  DOPLOC  system. 
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FIG.  5- LOCATION  OF  SATELLITE  RELATIVE  TO  CENTER  BEAM 


using  the  transmitter  at  Fort  Silli  Oklahoma,  and  the  receiver  at 
Forrest  City,  Arkansas.  The  transmitter  operated  continuously;  the 
receiver  was  turned  on  at  one  second  intervals  and  operated  until  1000 
cycles  of  Doppler  were  counted.  The  recorded  data  consisted  of  the  time 
duration  of  these  1000  cycle  Doppler  intervals,  together  with  the  starting 
time  of  each  Interval  (figures  6a  and  6B}.  The  average  frequency  shift 
over  each  1000  cycle  Doppler  interval  was  computed  from  the  data  by  the 
formula 


Afi 


I  ^ 


A-t- 


where  is  the  duration  of  the  Doppler  interval,  and  the  quantity  7000 

represents  an  empirical  constant  of  the  particular  transmitter-receiver 
system.  This  was  converted  to  the  range  rate  by 

pt  “  *2  ^^6 


and  was  assigned  to  be  the  value  of  p  at  the  midpoint  of  the  associated 
Doppler  interval  Atj  . 

The  preliminary  orbits  obtained  are  given  in  figure  7.  For  compari¬ 
son,  the  fined  solutions  obtained  by  BRL  and  NSSCC  (Space  Track)  are  also 
listed.*  BRL  uses  an  analog  procedure  to  determine  a  preliminary  orbit 
from  the  DOPLOC  data,  while  at  Space  Track  all  available  data  is 
compiled  over  a  period  of  time  to  compute  the  orbit. 


*Space  Track  results  were  furnished  to  the  author  by  BRL. 
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FIGURE  6  B 
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3.  ESTABLISHING  A  PRELIMINARY  ORBIT  FOR  ACTIVE  SATELLITES 


When  the  transmitter  is  aboard  the  satellite,  ^  2  O  and  the 


Doppler  equation  for  a  single  receiver  is 


(3.  1) 


Integration  of  this  equation  discloses  that  at  any  given  time  t  the  satellite 
lies  on  the  surface  of  a  sphere  of  radius  |0  (i)  ,  with  center  at  the 
receiver. 

Eliminating  the  subscript  R  for  convenience,  the  vector  form  of 
(3.  1)  is 

and  again  neglecting  the  earth's  rotation, 


A'f 

r 


(3.  2) 


where  V  is  the  satellite  velocity. 

Now  define  the  tinrje  of  closest  approach  as  the  time  at  which  ^  is 
minimum.  Then  the  following  conclusions  have  been  established  for  an 
active  satellite; 

At  the  time  of  closest  approach  t^  the  satellite  lies  on  the  surface  of 
the  sphere  of  radius  ,  with  center  at  the  receiver  R;  further,  the 

satellite  trajectory  is  tangent  to  this  sphere. 

If  we  now  assume  that  Tj*  2  p  is  constant  in  the  neighborhood  of 
closest  approach  t^  ,  then 


See  appendix. 
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and,  in  analogy  to  (2.  13), 


Thus,  in  the  case  of  the  active  satellite,  the  parameters  at  closest 
approach,  >  are  related  to  the  Doppler  data  by  (3.  1) 

and  (3.  3).  These  parameters  are  computed  in  the  same  way  as  were  the 
passive  satellite  parameters  /^/i  • 

The  procedure  for  determining  the  satellite  position  and  velocity 
components,  at  closest  approach  is  quite  analogous  to  that  for  the  passive 
satellite,  except  that  now  the  satellite  speed  at  closest  approach  is 
known  (  V  -  VI  )• 

If,  corresponding  to  the  base  line  TR  of  the  ellipse  of  figure  3,  we 
define  the  parameters and  for  the  circle  of  figure  8  with  respect 
to  the  latitude  of  the  receiver  R  ,  then,  as  before. 


<1 

<!><  ’'  =  TT-  J  yt  <  -X  <  TT^  , 


We  will  place  the  origin  of  the  topocentric  coordinate  system  at.  R  ,  with 
x-axis  vertical,  y-axis  due  North,  and  z-axis  due  West.  Then  the  first 
quadrant  of  the  circle  is  the  north-east  sector,  the  second  quadrant  the 
north-west  sector,  and  so  forth  (figure  8).  Thus,  the  relationships  that 
exist  between  the  computed  values  of  0  ,  ,  and  ,  and  the 

quadrants  in  which  the  satellite  might  lie  are  again  given,  by  figure  4. 


The  active  satellite,  at  altitude  x  »  will  lie  on  the  circle 

4  ,  _z  ^2 


Since  the  satellite  path  is  tangent  to  this  circle,  the  coordinates 
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FIG.8-0RIENTATI0N  OF  SPHERE  AND  ANTENNA  BEAMS 


o£  the  satellite  are  determined  from  the  expressions 


5/0*  o<  . 

As  before,  the  velocity  components  are  obtained  from 

K-O^  ^  ®  S'"?  cs<^  j? Ss  ^'■  CCS  oi  . 

The  following  conclusions  are  then  established  for  the  active 
satellite; 

Given  ^  t  >0.  altitude  X  and  G  ~  ar'C  nin 

/  •'  Vdii  ’ 


A  =  cot  28°  -  tan  4°^ 

and  denoting 


D, 

cos(eyLc)^ 

n 

r 

F. 

Ij 

Z\ 

ccs(&yLij^ 

sin  (Oyc) 

If  0<  yii. 

Cd.>c,  P^<c, 
^  >0 


,  there  are  two  possible  points  in  1st  quadrant  and  two 
possible  points  in  3rd  quadrant.  Since  '^/z  for 

all  these  points,  the  general  direction  of  motion  of  the 
satellite  is  known. 


Points  in  Ist  quad:  y  =  Dp  z  =  D^i  y  =  T  V  sin  (  Gyuf^  z  =  ±  a'"  cos  (  &  ) 

y  =  Fp  z  =  -F2;  y  =  i  'V  sin  (Gtyi'')>  z  =  ^  V  cos{  Gy-t) 

Points  in3rd  quad:  y  =  -Dp  z  =  -D2i  y  =  ^  T*  Bin{&  yu) ,  z  =  ±  a''cos(  6yc  ) 

y  =  -Fp  z  =  F2;  y  =  ±  sin(  ),  z  =  ±  V  coa{  Qyi^) 


* 

The  sign  of  y  is  to  be  determined  from  the  data  for  each  satellite. 
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\i  yx<  Q<  TL-yu.  •  there  is  one  possible  point  in^  each  o£  the  four  quadrants. 

(  X>,  >0^  P  >0  The  direction  of  motion  of  the  satellite  is  not  known.* 

FT  >'■■!)  R>o), 

Point  in  Ist  quad:  y  =  Fj,  z  =  -F2;  y  =  ±  V  sin  {  i  =  ±'xr  cob  ( 

Point  in  2nd  quad;  y  =  Dj,  z  =  D2;  y  =  *  V"  sin  (  ),  z  =  sp  V"  cos  (  ) 

Point  in  3rd  quad:  y  =  -Fj^,  z  =  F2;  y  =  ★  V  sin  (  Q  ■^),  i  =  ±  Vcob  {  B-pn- ) 
Point  in  4th  quad:  y  =  -Dj^,  z  =  "^2'  V  Bin(&^  ),  z  =T  V  cos  )  j 


If  K.-u.  11 

2  /  2 

fi>o) 


,  there  are  two  possible  points  in  the  second  quadrant 
and  two  possible  points  in  the  4th  quadrant.  Since 
ylz<0  for  all  these  points,  the  general  direction  Of 
motion  of  the  satellite  is  known. 


Points  in  2nd  quad;  y  =  D^,  z  =  D2;  y  =  d:Trsin(^^  ),  z  *  T  V  cos 

y  =  -Fj,  z  =  F2;  y  =  *  V  8in(  0^),  z  =  ;:t?  V  co8{  ^ 

Points  in  4th  quad:  y  =  -Dj^,  z  =  -D2;  y  =*V"  sin.  {  0 y*c)»  z  =  T  \rcoa{  ) 

y  =  Fj,  z  s  -F2!  y  =  ds  V  8in(  ,  i  =  ±  V^cos  (  &p^  ) , 


The  ambiguity  of  direction  in  the  second  case  (yu.<  &< 
removed  if  yU,  -  ^  ,  i.  e. ,  if  the  transmitter  antenna  is  rotated  45*^  from, 
the  latitude  of  the  receiver. 

Finally,  the  location  of  the  satellite  relative  to  the  center  beam  can 
be  used  as  a  criterion  for  selecting  permissible  points  (  y,  z  )  at  a  given 
altitude  x  ,  from  the  four  points  just  computed.  As  in  the  case  of  the 
passive  satellite,  this  criterion  is  given  by  inequality  (2.  35). 


s 

Although  the  sign  of  y  is  known  from  the  data,  the  sign  of  z  depends 
upon  whether  the  satellite  is  either  in  the  let  or  3rd  quadrants,  or  in  the 
2nd  or  4th  quadrants.  This  is  not  known  in  this  case. 


-  34  - 


4.  BIBLIOGRAPHY 


1.  R.  B>  Patton,  Jr. :  A  Method  of  Solution  for  the  Determination  of 
Satellite  Orbital  Parameters  from  DOPLOC  Measurements, 

BRL  Memo.  Report  #1237,  September  19S9. 

2.  W.  H.  Guier  and  G.  C.  Weiffenbach:  Theoretical  Analysis  of  Doppler 
Radio  Signals  from  Earth  Satellites,  John  Hopkins  University, 

Applied  Physics  Laboratory,  Bumblebee  Series  Report  No.  276, 

April  1958. 

3.  T.  Skinner:  Tracking  of  a  Moving  Transmitter  by  the  Doppler  Effect, 
AFCRC-TR-58.364,  (ASTU  Document  AD208699),  December  1958, 

4.  V.  A.  Kotelnlkov:  Determination  of  the  Elements  of  the  Orbit  of  a 
Satellite  Using  Doppler  Effect;  translated  from  Radiotekhnika 
ElektHckestvo,  (7),  1958;  Reproduced  by  the  National  Academy  of 
Sciences. 

5.  E.  G.  C.  Burt:  The  Computation  of  Orbit  Parameters  from 
Interferometer  and  Doppler  Data,  Proceedings  of  the  Royal  Society, 

28  October  1958. 

6.  R.  R.  Newton:  Tracking  Objects  within  the  Solar  System  Using  Only 
Doppler  Measurements,  Proceedings  10th  International  Astronautical 
Congress,  Volume  1,  1959* 


-  35  - 


5.  APPENDIX 


Derivation  of  the  Doppler  Equation  for  a  Passive  Satellite. 


We  will  assume  that  radio  waves  travel  with  constant  velocity  of  light 
through  all  disturbances,  and  that  the  motion  of  the  satellite  is  slow  com¬ 
pared  to  the  speed  of  light. 

We  consider  a  stationary  transmitter  T  on  the  surface  of  the  earth 
at  sonrie  positive  distance  from  a  stationary  receiver  R  on  the  surface 
of  the  earth.  A  signal  emitted  at  T  at  time  t,  travelling  at  the  speed 
of  light  c  ,  arrives  at  the  moving  satellite  S  at  time 

•c=  t  , 

where  ^  is  the  magnitude  of  the  position  vector  of  the  satellite  S 
relative  to  the  transmitter  T  .  Denote  by  an  active  transmission 
time  Interval  and  by  ^  the  corresponding  time  interval  during  which 
the  satellite  is  Illuminated.  Then 

p,  •  Pt  ^  g  ^ —  • 

If  the  relative  range  of  the  satellite  Is  assumed  to  be  a  linear  function 
of  time  during  this  interval,  we  obtain,  dividing  by  , 

I  =  ^ 

U 

which  yields 

IR  \\ 

'  -7-  >  (5.  1) 

where  and  'f^  are  the  reflected  (satellite)  and  transmitted  frequencies, 
respectively. 

In  a  similar  manner,  the  relationship  between  the  satellite  illumination 
interval  and  the  interval  during  which  the  reflected  signal  is  received  at  R 
is  expressed  in  terms  of  the  satellite  frequency  %  ^be  receiver 
frequency  by  the  equation 

i  =,  I  -  i  ,  (5.  i) 
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where  higher  order  terms  in  fg/c  have  been  neglected.  From  (5.  1)  and 
(5.  2)  we  obtain  the  following  expression  for  the  Doppler  frequency  shift 
of  a  passive  satellite: 


and,  again  neglecting  second  order  terms  in  ^/c  ,  this  simplifies  to  the 
Doppler  equation  for  a  passive  satellite, 


(5.3) 


If  the  transmitter  is  aboard  the  satellite,  ^^'5  p  and  equation  (5.3) 
reduces  to  the  Doppler  equation  for  an  active  satellite, 
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